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Abstract 

A parameter- free geometric model for nuclear ab- 
sorption is derived herein from microscopic theory. 
The expression for the absorption cross section in 
the eikonal approximation, taken in integral form, is 
separated into a geometric contribution that is de- 
scribed by an energy- dependent effective radius and 
two surface terms that cancel in an asymptotic series 
expansion. For collisions of light nuclei, an expression 
for the effective radius is derived from harmonic os- 
cillator nuclear density functions. A direct extension 
to heavy nuclei with Woods-Saxon densities is made 
by identifying the equivalent half-density radius for 
the harmonic oscillator functions. Coulomb correc- 
tions are incorporated, and a simplified geometric 
form of the Bradt-Peters type is obtained. Results 
spanning the energy range from 1 MeV/nucleon to 
1 GeV/nucleon are presented. Good agreement with 
experimental results is obtained. 

Introduction 

The nuclear absorption cross section is a mea- 
sure of the sum of all the reaction processes induced 
during nuclear collisions. This metric is applied 
widely in such diverse fields as fundamental nu- 
clear physics, the study of radiation effects on living 
cells, and radiation-shielding design for future space- 
exploration vehicles (ref. 1). 

In recent years, considerable advancements in 
both the theory (refs. 2—19) and experimental meth- 
ods (refs. 20-26) of this fundamental observable have 
been made. Most significantly, the energy varia- 
tions of the absorption cross sections for all nuclei in 
the intermediate-energy range from 10 MeV/nucleon 
to 1 GeV/nucleon are qualitatively similar and 
can be related to the underlying two-nucleon inter- 
actions. At energies above 1 GeV/nucleon, the ab- 
sorption cross section is essentially constant and ap- 
proximately equals the total geometric cross-sectional 
area of the projectile-target system; this consistency 
allows an extension of the geometric model and its 
parameterized forms (refs. 13-17, 25, and 26) to the 
entire intermediate- and low-energy range. The geo- 
metric formula provides a quick method for the cal- 
culation of the absorption cross section and is highly 
desirable for inclusion in radiation-shielding trans- 
port codes. 

Microscopic formulations are based on the quan- 
tum collision theory of many-body composite sys- 
tems. The solution of this complex problem can be 
approached with the methods of Glauber (ref. 2), 
Watson (ref. 3), Feshbach (ref. 3), and the KMT ap- 
proach (ref. 4). Wilson (ref. 5) extended the Watson 


formalism of nucleon-nucleus collisions to nucleus- 
nucleus collisions, from which the cross sections in 
the eikonal approximation can be obtained. These 
solutions, which are generally solved by numerical 
methods (refs. 5 and 6), are given in integral form. 

Closed-form solutions for the microscopic theory 
have been obtained for the simplest case of Gaussian 
nuclear density distribution (refs. 7 and 8). The 
solution applicability to the more realistic density 
distribution of the Woods-Saxon form is based on 
matching equivalent Gaussian forms to the surface 
densities of heavier nuclei. This approach is suc- 
cessful because the energy dependence for the trans- 
parency is largely a function of the low-density sur- 
face region with a saturated, absorptive inner core 
that remains essentially energy independent. 

The purpose of this paper is to derive a geomet- 
ric formula, valid for both light and heavy nuclei, for 
the nuclear absorption cross section from the micro- 
scopic formulation. This effort will be accomplished 
by separating the integral expression for the absorp- 
tion cross section in the eikonal approximation into 
an effective geometric radius and two surface com- 
ponents that approximately cancel in the asymptotic 
expansion. 

The concept of the effective radius in microscopic 
formulations is not new (refs. 8 and 9). However, in 
all previous works the contributions beyond a cutoff 
radius are negligible, whether the effective radius 
is defined at half-absorption density (ref. 8) or at 
a distance of one mean free path from the surface 
(ref. 9). In this work, however, the effective radius at 
half-absorption density has a precise definition such 
that the transparency within this radius cancels the 
absorptive contribution from the outer region. 

A summary of the format for this paper is pre- 
sented. In “Theory,” the results from the micro- 
scopic theory are reviewed and developed into a geo- 
metric formalism. In “Application to Light Nuclei,” 
a method is illustrated for harmonic oscillator mat- 
ter densities that are applicable to light nuclei. Ex- 
pressions for the eikonal phase shift and the total 
absorption cross section are derived for the nucleon- 
nucleus as well as for the nucleus-nucleus inter- 
actions. An exact formula for the effective radius 
at half absorption is derived, and its first approxima- 
tion is used to evaluate explicitly the geometric term 
and the surface contributions. Results comparing ex- 
periments for energy values from 50 MeV/nucleon 
to 1 GeV/nucleon are presented. Also, a simpli- 
fied version of the expression for the effective radius, 
which is used in the rest of the sections, is given. 
In “Extension to Woods-Saxon Densities,” a direct 



extension of the expression mentioned above is made 
for the heavier nuclei described by Woods-Saxon den- 
sity distributions. In “Coulomb Interactions,” the 
effects of Coulomb deflections are incorporated, and 
results for light- and heavy-ion projectile-target sys- 
tems spanning the energy range from 1 MeV/nucleon 
to 1 GeV/nucleon are presented. In “Bradt-Peters 
Form,” the expressions are further simplified and the 
results compared with those of the previous section. 
In “Nuclear-Medium Effects,” comments are made 
on nuclear-medium effects. “Concluding Remarks” 
summarizes the suggestions for improvement submit- 
ted herein. 

Theory 

Microscopic Theory 

The nuclear absorption cross section in the eikonal 
approximation is given by 

roc 

Gu bs = 27 r / {1 - exp [—21m X(6, e)]}b db (1) 

JO 

where Im denotes the imaginary part of the complex 
phase shift X(6, e), which is a function of the impact 
parameter b and energy e, and can be derived from 
a multiple-scattering series expansion for the optical 
potential (refs, 5 and 6). Within the Watson form 
of impulse approximation, when the propagator G is 
replaced by the free propagator G 0 (ref. 5), the phase 
shift is given by 

X{b ’ e) = k(Apt T A T ) I dz I 


the real to imaginary components of the forward- 
scattering amplitudes. 

Geometric Model 

A geometric formula for the absorption cross sec- 
tion of high-energy nuclear collisions is given by the 
Bradt-Peters form (ref. 2), 

^abs = * (Rp + r T “ W 2 ( 4 ) 

where Rp and Rp are the radii of the projectile and 
target nuclei and 8 0 is a constant. If 8 0 is small com- 
pared with the radii, the formula shown above implies 
an annular region of transparency for peripheral colli- 
sions. To show that this form may indeed be derived 
from microscopic theory, the integral given in equa- 
tion (1) is separated at an impact parameter R e as 
follows: 

fRe 

<7 a bs = 27r / {1 - exp [-2ImX(6, e)]} b db 

Jo 

roc 

+ 27 r / {1 — exp [— 2ImX(6,e)]} b db 

J Re 

= nRl - Si + S 2 (5) 

such that <5i and S 2 are the small surface contribu- 
tions given by 

fRe 

6 \ = 2tt I exp [-2ImX(6, e)]b db (6) 
Jo 

and 


X JfPyPpib + z + y + $)t(e,y) (2) 6 2 = 2x / {1 - exp [-2ImX(f>, e)]} 6 db (7) 

J Rt 


where k is the wave number of the projectile, Ap 
and Ap are the mass numbers, and pp and pp are 
the projectile and matter densities; £(e, y) signifies 
the two-body transition amplitude for free nucleons 
as follows: 

t(e, y) = - (— ) ^ a(e) [a(e) + i] [ 2n B(e)]~ 3 / 2 
\m / 


x exp 


2 B(e) 


(3) 


Here, e is the two-nucleon kinetic energy in the 
center-of-mass (cm) frame of reference, m is the 
nucleon rest mass, a(e) is the isospin-averaged total 
nucleon- nucleon cross section, B(e) is the averaged 
slope parameter, and a(e) is the averaged ratio of 


Appendix A includes the asymptotic series 
expansions, 


as 2nR e 


exp[— Xi(fle,e)] 

dX l {b,e)/db\ Rt 


( 8 ) 


and 


h » ■ ( 1 ~r’i~vS ' e)l}2 “Pl*> <*. «» 

0X1(0, e)/db\R e 


where 


Xi (6, e) = 2ImX(6, e) 


(9) 


( 10 ) 


Note that 62 can also be expanded into a conver- 
gent Taylor series, provided that (6, e) < 1 beyond 



R e and that, in the first approximation, 62 is given 



From equations (8) and (9), 

62 — 6\ = 2nR € 6(R e} e) (12) 


where 


6(R e ,e) 


Hence, 


exp[Xi(fl e ,e)] , 
dX l (b,e)/db\ Rt { 


2exp[-Xi(P e ,e)]} 

(13) 


o'abs = ’r-Re ~ 2TrR e 8(R e ,e) 
« 7T [Re - 6{Re,e)} 2 


(14) 


If the integration limit R e is the sum of the projec- 
tile and target radii, the required geometric formula 
shown in equation (4) is obtained. However, for con- 
venience, R e can be chosen such that <5 is identically 
zero. When this condition is valid, equation (13) 
shows 

exp[— Xi(i? e , e)] = 0.5 


or 

Xi(flc,c) = In 2 (15) 

That is, Re is chosen at half transparency or, equiv- 
alently, at half absorption. (See fig. 1.) The trans- 
parency function T is 

T = exp[-Xi(&, e)] (16) 

and the absorption function A = 1 — T. With <5 = 0, 

^abs = fRe ( 17 ) 

where the effective radius R e is now a function of 
energy. This equation gives the desired geometric 
formula for the nuclear absorption; R e is derived 
from equation (15). However, a more appropriate 
formula for Re can be derived in terms of the classical 
geometric radius R 9 y which is defined as oo* 

Note that at the high-energy limit, as e — ► oo, the 

wavelength of the projectile nucleons is sufficiently 
small compared with the nuclear dimensions, and the 
effective radius assumes the classical geometric radius 
Rg, given by 


Rg = To ( A J / 3 + 4 /3 ) (18) 

where r 0 is a variable determined by the condition 
Xi(R g , oo) = In 2 (19) 


such that 5 = 0 is satisfied. The effective radius R e 
at an arbitrary energy e can then be related to the 
geometric radius R g by equating 


X x ( R e , e) = X x ( R g , oo) = In 2 (20) 

This relation is useful only if an algebraic expres- 
sion for X(6, e) is known. For Gaussian and harmonic 
oscillator density functions, X(6, e) is readily obtained 
by the methods in “Application to Light Nuclei.” 

Application to Light Nuclei 

In this section, the phase-shift function X(6, e) is 
derived from the harmonic oscillator density func- 
tion, which is applicable to light nuclei for A < 16. 
Also, an exact expression for the effective radius is 
obtained. The nucleon- nucleus and nucleus-nucleus 
collisions are treated separately. 

For the nucleon-nucleus collision, the nucleon 
density is simply 


p(x - x f ) = 6 {x - x') (21) 

because the nucleon wave function has been effec- 
tively incorporated into the two-body transition am- 
plitude £(e,y), and the nucleus charge density is 
given by 

Pc(x) = Po ^ + 7^2 ) ex P ^2 j ( 22 ) 

where 7 and a fm (1 fermi —lx 10~ 15 m) are 
constants taken from reference 27. The expression for 
X(6, e) is derived in appendix B and has the following 
form for the nucleon- nucleus system: 

Xi(6, e) = /i(e) (ai + a 2 & 2 ) exp[-P(e)6 2 ] (23a) 

and for the nucleus-nucleus system: 

Xi(6, e) = /i(e) (ai + a 2 & 2 + a^b 4 ^ exp[-P(e)6 2 ] 

(23b) 

where explicit expressions for the functions f\ (e) and 
P(e) and for the constants ax, 0-2, and 03 in each case 
are given in appendix B. 

Figure 2 displays the absorption function 

A — 1 — exp[~Xi (by e)] (24) 

as a function of the impact parameter 6, for 4 He- 4 He 
and 16 0- 16 0 systems, at an energy of 2 GeV/nucleon, 
which is taken as the geometric limit. The position of 
the effective radius Re that is defined at T = A = 0.5 
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is also indicated. Note that even within the geometric 
limit, considerable transparency for peripheral colli- 
sions occurs. 

The effective radius R e for three different energy 
values is shown in figure 3, where the absorption 
is plotted at 50 MeV/nucleom 250 MeV/nucleon, 
and 1 GeV/nucleon for the ™0- l6 0 system. As 
shown in the figure, Re passes through a minimum 
and then increases as the energy decreases. The 
strong absorption radius R s is reached in the low- 
energy limit. The underlying basis for this energy 
dependence can be traced to the nucleon-nucleon 
total cross section, which exhibits a similar trend. An 
expression for the effective radius R e in terms of its 
geometric value Rg is obtained from equations (23) 
and the condition 

For nucleus-nucleus collisions, 


Table I. Collision Values for r G 


Atomic mass 
number 

r Oi 

Fm 

Nucleon-nucleus 

Nucleus-nucleus 

2 

0.50 

0.47 

4 

1.02 

0.85 

7 

1.14 

1.13 

9 

1.19 

1.17 

10 

1.19 

1.16 

11 

1.18 

1.14 

12 

1.19 

1.13 

14 

1.20 

1.16 

16 

1.23 

1.21 

>16 

1.25 

1.25 


in figure 4 and listed in table I. From equation (20), 
this determination is clearly sufficient to ensure that 
the condition 6 = 0 is also satisfied by the value of 
the effective radius given by equation (25). 


Rl = P(00) * 2 


P(e)^ + P(e) \ ln 


' fl(e) ' 

./ l(oo). 


f Ol +02 Pe +03^ \ 1 

\ai + 02-R^ + a 3 R* J j 


(25) 


with a similar expression for the nucleon-nucleus case 
for a 3 = 0. In the equation above, the geometric 
value R g is given by 


Rg = 

(26a) 

for the nucleon-nucleus system and 


R g = r T Af + rpAp 3 

(26b) 


for the nucleus-nucleus system. The value r is deter- 
mined from the condition 


X(Rg } 2A GeV) = In 2 

which, for the nucleus- nucleus system, is given by 


For the heavier nuc lei within the harmonic oscil- 
lator category, r 0 « 1.2 fm; this length approximates 
the generally accepted value that determines the nu- 
clear radius. Clearly, the values of r Q for the lightest 
nuclei such as deuterium and helium are considerably 
smaller than those for the heavier nuclei; this differ- 
ence can be explained by examining the convolution 
integral for the phase shift. In the case of Gaussian 
density distributions, the resultant width increases 

by the factor (ap + j , whereas the total geo- 
metric width is given by the sum ap + ap, as shown 
in equation (26), where ap and ap are the widths 
of the target and projectile, respectively. Hence, for 
the lightest nuclei, which have near-Gaussian distri- 
butions, r 0 must be lowered to account for this dis- 
crepancy. The convolution widths for the Woods- 
Saxon forms, however, are expected to be close to 
the geometric value. 

The advantage of writing R e in terms of Rg } as 
given in equation (25), is that in the first approxima- 
tion the second logarithmic term may be neglected 
when R e & R gi with the result 


In/, ( 2 A GeV) + In (01 + a 2 R 2 g + a 3 flj) 

- P(2A GeV)P 2 = ln(ln 2) (27) 

Substituting for R g from equation (26), equa- 
tion (27) is solved iteratively with an initial value 
of r a = 1.0 fm. For the nucleus-nucleus system, iden- 
tical taxget-projectile systems are considered with 
r T = rp = r Q . The resulting values for r 0 are plotted 


Ri 


P(2i4 GeV) 2 , 1 

P(e) 9 P(e) 


In 


Me) 


fl(2A GeV). 


(28) 


Although this approximation offsets the condition 
6 = 0, the integration limit in equation (5) is suffi- 
ciently defined to evaluate 6] and S 2 explicitly; the 
absorption cross section is given accurately by 


°abs = nRe - S\{Re, e) + S 2 (R e , e ) (29) 
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Figures 5-9 display the absorption cross section 
using equation (29) for light-ion projectile-target sys- 
tems that agree well with experimental data. These 
results are also in close agreement with the numerical 
evaluation of the integral given by equation (1). (See 
refs. 5 and 6.) 

Hereafter, the expression for R e from equa- 
tion (28) will be used with the assumption that 
(#2 — ^i) is negligible. In equation (28), the expres- 
sion P(e) for the nucleon-nucleus collision is equiva- 
lent to equation (BIO), 


two-parameter Woods-Saxon charge density is given 
by 

Pc(r) = y - — 

i + exp(^) 

where R is the radius at half density, c is related 
to the nuclear skin thickness t by the expression 
c = i/4.4, and p 0 is a normalization constant. From 
reference 27, the parameter c is a constant approx- 
imately equal to 0.55 fm for all nuclei. Hence, we 
assume that only one parameter R describes the 
Woods-Saxon form. 


P(e) = 


a T ~ \ T p + 2B ( e ) 


(30a) 


The equivalent half-density distance R for the 
harmonic oscillator function satisfies 


and 5(e) for nucleus-nucleus collision is equivalent to 
equation (B17), 



P(e) = 


a T + ap — -r 2 + 25(e) 


-l 


(30b) 


which can be rewritten as 

k = In 2 + ln(l + 7 k) 


(34) 


where r p = 0.86 fm is the radius of the proton. 
Because the slope parameter 5(e) is approximately 
a constant equal to 0.4 fm 2 within the energy range 
considered, P(2A GeV) £3 5(e); equation (28), with 
R g given by equation (26), then simplifies to 

5 2 = + ( 0 %? + 0.35^ g(e) (31a) 

for the nucleon- nucleus collision, and 

2 

5 2 = (VpAp 3 + + ( a \ + a 2 p ) g(e) (31b) 


where 

d2 

k = (35) 

Equation (34) is solved iteratively for k with 
an initial value of k = In 2. The parameter k is 
plotted in figure 10 as a function of atomic mass 
number. For the heavier nuclei with the harmonic 
oscillator density form, k « 2 or 5 « y/2 a. Also, 
for these nuclei an equally valid description would be 
the Woods-Saxon form with the half-density radius 
given by R = y/2a. This equivalence suggests that 
the substitution of 


for the nucleus- nucleus collision, where 
g(e) = In 


a(2A GeV) 


Me) 

./1 (2 A GeV). 
r 9(e) 1 


R 



into equations (31) can extend its range of applicabil- 
ity to the Woods-Saxon densities of heavy ions with 
A > 16. This substitution yields for the nucleon- 
nucleus system, 


The terms Ap and Aj are the charge density width 
parameters for the harmonic oscillator function given 
in reference 27, r P and r p are the respective values 
for r 0 given in table I, and the parameterized form 
for a(e) is given in appendix C. 

Extension to Woods-Saxon Densities 

The Woods-Saxon density does not yield a simple 
expression for the phase shift X(6, e). Therefore, a 
simple modification of equations (31) is given. The 


R\ = (r T A l f) 2 + + 0.35^ 5(e) (36a) 

and for the nucleus-nucleus system, 

R\ = (r P A),/ 3 + r T Ay 3 ) 2 + \(R 2 t + R 2 p ) 9(e) 

(36b) 

where Rp and Rp are the half-density values for the 
Woods-Saxon forms given in reference 27. 
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Coulomb Interactions 

The derivation of the eikonal approximation as- 
sumes straight-line trajectories for the projectile and 
an impact parameter at the collision site that is 
equal to its asymptotic value. For heavy ions at 
low energy, the Coulomb force deflects the straight- 
line trajectories significantly, thereby decreasing the 
cross-sectional area of incident particles that actually 
undergo nuclear forces. The eikonal approximation 
may still be used (refs. 8, 10, and 18) if the effective 
impact parameter at the collision site is the distance 
of closest approach d (refs. 8 and 18), where b and d 
are related by 

*-*(>- 7 ) <37) 

with 

_ Z T Z P e l 
a 

^cm 

where Z t and Zp are the charges of the target and 
projectile, respectively, e Q is the electron charge, and 
E cm is the total center-of-mass energy. 

Hence, the asymptotic impact parameter b ei cor- 
responding to the radial distance Re , is given by 

*-* 0 -!) (38) 

where b e now defines the integration limit for equa- 
tion (5). Therefore, the Coulomb- modified geometric 
formula is readily obtained as 

<M» - (1 - j£) (39) 

In figures 11-21, the absorption cross section 
for several light- and heavy-ion projectile-target nu- 
clei is plotted with and without the Coulomb cor- 
rections for energy values from 1 MeV/nucleon to 
2 GeV/nucleon. The results with Coulomb correc- 
tions agree well with experiment. The following val- 
ues for r 0 (see table I) were used for nucleus-nucleus 
collisions: 

A = 4; r 0 — 0.56 fm 
7 < A < 16; r Q = 1.20 fm 

A > 16; r Q — 1.25 fm 

Figures 22-26 show the cross sections for spe- 
cific projectiles on various targets with and without 
Coulomb interaction. The linear relationship with 
atomic number is clearly evident. As expected, the 


Coulomb corrections are more significant for heavy 
ions. 

Bradt-Peters Form 

In the present formalism, the geometric limit for 
the nuclear absorption cross section is given by 

<7 a bs = * 1 % ( 4 J/ 3 + 4 /3 ) 

compared with the Bradt-Peters form at high energy, 

<7 a bs = 7rr o (^p 3 + - So) 

The difference comes from the choice of r Q . Here, 
r 0 satisfies the geometric criteria that are valid at 
high energies and small wavelengths such that 

R g = To ( A }/ 3 4- aJ/ 3 ) (T = 0.5; c -» 00 ) 

An equally valid description is to define the strong 
absorption radius R s at low energy as the reference 
value such that 

R s = r s (A l J, Z + 4 /3 ) (T = 0.5; e -> 0) 

The strong absorption parameter r s is greater 
than r 0 because the absorptive content of the surface 
region increases at lower energies due to s-wave scat- 
tering. This low-energy description results in 6 0 / 0 
at high energy, a state that conforms to the Bradt- 
Peters formula. 

However, since the eikonal approximation is valid 
only at high energies, the low-energy description 
cannot be implemented. Also, the corrections that 
show the Coulomb and Pauli correlation, significant 
contributors at low energies, were not incorporated 
into the phase-shift function. Hence, the choice of 
the geometric limit as the reference point is justified. 

The more general equations (36) for the effective 
radius R e at intermediate energy will now be cast in 
the Bradt-Peters form. The simplest form is obtained 
for heavy ions when rp = rp — r 0 . Using the 
parameterized form for half-density charge radius R , 
given by 

R — c\A 1 ^ - C2 

with 

ci = 1.18 i 

(40) 

C2 = 0.48 J 
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in equations (36) and simplifying, the Bradt-Peters 
form is obtained for the nucleon-nucleus system 


R 2 e=r 2 (e)[ 


.1/3 


60(e) 


where 


r 2 (e) = r 2 + 


c? 9(e) 


(41a) 


6M “ 

and for the nucleus-nucleus system as 

R 2 e = r 2 (e) [>lp /3 + Ap /3 - 6 0 (A P ,A T ,e)] 2 (41b) 

with r as shown before and 


in the intermediate-energy region. At high energy, 
many scattering states are available; hence, the ef- 
fects of the occupied states are comparatively negli- 
gible. At low energy, peripheral collisions dominate 
when the nuclear overlap densities are low and when 
the number of occupied states is small. Peripheral 
collisions dominate at low energies for the following 
reasons: the nuclear absorptive power increases due 
to s-wave scattering between nucleons, and Coulomb 
forces deflect the trajectories into the peripheral re- 
gions. Thus, reasonably good results are obtained 
with the Coulomb corrections alone. The nuclear- 
medium effects should actually increase the mean 
free paths (ref. 19); hence, more nuclear transparency 
results. For uncharged projectiles such as the neu- 
tron, these effects can be considerable, as shown in 
figure 14(b). 


<5o(Ap, Ap,e) « 


9(e) 

2 r 2 (e) 


^cic 2 +c\ 


A l/3 A 1/3 \ 

aP + a'*)' 


where r 0 = 1.25 fm. 

In deriving the formula above, 8 0 is assumed to be 
1/3 1/3 

much smaller than Ap and Aj, . Note that r is now 
energy dependent, which supports the observation 
made by Kox et al. (ref. 25). 

After the above equation is substituted into equa- 
tion (39) to account for Coulomb correction, the re- 
sults are compared with those of the previous section 
in figures 27-30. The agreement is good for the heav- 
ier as well as for the proton projectiles. The discrep- 
ancy for the alpha projectiles in the nucleus-nucleus 
case results from the assumption rp — rp = 1.25 fm 
for all nuclei in the Bradt-Peters form. 


Nuclear-Medium Effects 

In the microscopic formulation, the two-nucleon 
interaction is defined between free nucleons. This 
definition excludes nuclear-medium effects from 
Fermi motion and Pauli blocking. The exclusion 
of these effects can be justified in both the high- 
and low-energy (below 20A MeV) regimes but not 


Concluding Remarks 

A parameter- free geometric model for nuclear ab- 
sorption was derived herein from microscopic theory. 
The expression for the absorption cross section in the 
eikonal approximation, taken in integral form, was 
separated into a geometric contribution that was de- 
scribed by an energy-dependent effective radius and 
two surface terms that cancelled in an asymptotic 
series expansion. 

For collisions of light nuclei, an expression for the 
effective radius was derived from harmonic oscilla- 
tor nuclear density functions. A direct extension to 
heavy nuclei with Woods-Saxon densities was made 
by identifying the equivalent half-density radius for 
the harmonic oscillator functions. Coulomb correc- 
tions were incorporated, and a simplified geometric 
form of the Bradt-Peters type was obtained. Re- 
sults spanning the energy range from 1 MeV / nucleon 
to 1 GeV/nucleon were presented. Good agreement 
with experimental results was obtained. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
April 23, 1993 


7 


Appendix A 

Surface Contributions to Total Absorption Cross Section 

The surface contributions to the total absorption cross section, given by <5i and 62 , are evaluated using the 
Laplace method for asymptotic series (ref. 27) 


fRe 

6\ =2n / exp[— X\(b, e)]b db 
Jo 

= 2 -K J exp|- jxi(i? e ,e) + — {*> - R e ) + . ..j j 6 db 

- p exp[-Xi(fle,e)] 

dXi (b,e)/db\ Re 


(Al) 


Similarly, 


/‘OO 

62 = 2 n / {1 - exp[-Xi(6,e)]}fe db 

jRe 

rOO 

= 27T / exp(ln {1 — exp[— X\(b, e)]} )b 

J Re 

= 2tt / c 
JRe 


db 


ln{l - exp[— Xi(/Z et e)]} + — (ln{l - exp[-Xi(6, e)]}) (b - Re) + .. J b db 


2nR e 


{1 - exp[-Xi(fl e ,e)]} 2 
exp[X i (R e , e)] [3Xi (6, e) / db] | ^ 


(A2) 
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Appendix B 

Derivation for Phase Shift X(6, e) 

Expressions for the phase shift X(b,e) are derived from the harmonic oscillator density functions, 
charge density is 


Pc{ x ) — Po 



The 


(Bl) 


where 7 and a fm are taken from reference 27 and p a is a normalization factor. 
The matter density is derived from the charge density (ref. 29) as 


Pm{x) = Po{o! + b'x 2 ) exp(-px 2 ) 


(B2) 


where 


a '- 8?[ 1 + 


3 3 a 2 \ 

2 7 "8Vj 


b' = 7 


128s 7 


(B3) 

(B4) 


where r p — 0.86 fm is the radius of the proton, 
separately. 



P ~ 4s 2 

Nucleon-nucleus and nucleus-nucleus collisions are considered 


N ucleon-N ucleus 

The matter density function for the nucleon is given by 

p p (x - x') = 6(x - x') 

and for the target nucleus by 

Pr(x) = Pro( a T + b' T x 2 ) exp {-px 2 ) 

The phase shift is given by equation (2), 

*( b , e ) = ]7 (ap+ T a t ) j dz j J d3 ypp( b + z + y + y) 

Substituting for pp, pj, and t(e,y) from equation (3), the integral is easily performed, obtaining 

Xi(6,e) = 2ImX(6, e) = /i(e)/ 2 (6,e) 

where 


/l(e) = A T p TO a(e) 


7T 


1/2 


[2B(e)(p + p')}- 3/2 


P(e) J 

/ 2 (6, e) = (ai + a 2 b 2 ) exp[-P(e)6 2 ] 


(B5) 

(B6) 


(B7) 

(B8) 

(B9) 
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i III l|lll I IN III II II II. 


with 


, „ _ b r a,2 

cl\ = a + 1.5 —j 4- 


a 2 = 


p + pf 2 P(e) 

f 2 

p b> 


P = 


P = 


(p + p') 2 
1 


2 B{e) 
1 

4sj, 


P(e) 4 4 + 25(e) 

Also an expression for <5 2 in the Taylor expansion is derived as 

1 3/2 


62 = 2 


7T 


LP(e)J 


/i(c) 


fll + L5 P^) +a2 ^ 


exp 




N ucleus-N ucleus 

The matter density function for the projectile is 


and for the target is 

which gives 

where 

and 

with 


pp(x) = ppo (a!p + b' p x 2 ) exp(-px 2 ) 

pp(x) = Pro( a T + b T x 2 ) exp (-qx 2 ) 

Xj ( 6 , e) = /i(e)/ 2 ( 6 ,e) 

/l(e) = ApAp d{e)pTOPPO {27r5(e)[p + pVJHfl + g 7 (e)]} 3 ^ 2 
/a( 6 , e) = (ai + a 2 6 2 + 03 6 4 ) exp[-P(e) 6 2 ] 


p - 


4sp 


4s 2 


9' = 


1-1 


25(e) 

p'{e) = ^4sp + 25(e)j 
P(e) = [4 ^sp + 4) + 25(e)j J 


(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 
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dl = 


a P a T + ^^7 ( a T b P + a P 6 r) + - 4 - ( p + ^)2 

//./ P' 2 P 2 . / - / « PP' 1 ./ w 3 P 2 +P' 2 

= <>T b P + a P b TJ^Yjl ~ b Tbp2j~^ + b T b p 2 {p - p , )3 


_ a ./ P 2 P' 2 

a 3 — °T°P] 


(p + p 7 ) 4 


where a! p and 6 p are the constants that relate to the projectile as given by equations (B3) and (B4) and 

, 3 _ Zl T a 2 T \ 7 r 4 15 

\ + 2 7r 8s p y 16s p ( q + q ')_ 


„3 

T 8s| 


. _ _a|_ 7 t 4 

r 8s| 164 


'2 


(<7 + <?')' 


ofp — dp 4“ bp 


1.5 

q + q' 


relate to the targets. 

An expression for 62 in the Taylor series expansion is derived 


as 


62 


Me) 


/ 7T f 

1.5 d2 ‘ 


lm\ ai+a 2 

[p (e) + 

+ 03 


15 


ml 

4 P(e) 2 + P(e) 


+ 4Ri 


jexpj- 


P(e)R; 


The following integral relations were used. 

1 = /” expI “ ,r2 ■ p{r ~ y,2] * = (jT^) ' Btp (~;Ti ) yl 


(bis) 


fRicn 


r\ 

J r 2 exp[— gr 2 - p(r - y) 2 ] dr = — —I (B20) 

J{ r - y) 2 exp[— gr 2 - p(r - y) 2 ] dr = - (B21) 

J r 2 {r-y) 2 exp[-qr 2 -p(r-y) 2 ]dr= (~^) 7 ( B22 ) 

An alternate partial evaluation is given in reference 30. 
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Appendix C 

Parameterization for Experimental Two-Nucleon Total Cross Section 

For completeness, the isospin-averaged nucleon-nucleon cross section taken from reference 1 is 
parameterization for the experimental two- nucleon total cross section in millibarns (1 barn = 1 x 10 
given for the proton-proton interaction by 

crpp(e) = (l + h - ) {40 + 109 cos (0.199\/e) exp [-0.451 (e - 25) 0 ' 258 ] } 

for e > 25 MeV and for the lower energies by 

crpp(e) = exp |6.51 exp |- (^)° J | 

For proton-neutron interaction at e > 0.1 MeV, 

a np (e) = 38 + 12 500exp[— 1.187(e - O.l) 0 35 ] 

and at lower energies 

a np (e) = 26000exp[-(^)°' 3 ] 

The isospin average is calculated as follows: 

<r(e) = [(NpNp 4 ZpZp) Gpp 4 ( NpZp 4 ZpNp) a np ] JApAp 

where A } N } and Z are the atomic, neutron, and proton numbers, respectively. 

The nucleon- nucleon slope parameter is parameterized (ref. 31) as 

B(e) = 10 4 0.5 In ( — 

\s 0 

where s f is the square of the nucleon-nucleon center-of-mass energy and s 0 = (1 GeV/c)~ 2 . 


given. A 
“ 28 m 2 ) is 

(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 


12 


References 

1. Wilson, John W.; Townsend, Lawrence W.; Schimmerling, 
Walter; Khandelwal, Govind S.; Khan, Ferdous; Nealy, 
John E.; Cucinotta, Francis A.; Simonsen, Lisa C.; Shinn, 
Judy L.; and Nor bury, John W.: Transport Methods and 
Interactions for Space Radiations . NASA RP-1257, 1991. 

2. Glauber, R. J.: High-Energy Collision Theory. Lectures 
in Theoretical Physics , Volume J, Wesley E. Brit tin and 
Lita G. Dunham, eds., Interscience Publ., Inc., 1959, 
pp. 315-414. 

3. Joachain, Charles J.: Quantum Collision Theory , Third 
ed. Elsevier Science Publ. Co. Inc., 1983. 

4. Kerman, A. K.; McManus, H.; and Thaler, R. M.: The 
Scattering of Fast Nucleons From Nuclei. Ann. Phys . 
(N.Y.), vol. 8, no. 4, Dec. 1959, pp. 551-635. 

5. Wilson, John W.: Composite Particle Reaction Theory. 
Ph.D. Diss., College of William and Mary in Virginia, 
June 1975. 

6. Wilson, J. W.; and Townsend, L. W.: An Optical Model 
for Composite Nuclear Scattering. Canadian J Phys., 
vol. 59, no. 11, 1981, pp. 1569-1576. 

7. Karol, Paul J.: Nucleus-Nucleus Reaction Cross Sections 
at High Energies: Soft- Spheres Model. Phys. Review C , 
vol. 11, no. 4, Apr. 1975, pp. 1203-1209. 

8. Charagi, S. K.; and Gupta, S. K.: Coulomb-Modified 
Glauber Model Description of Heavy-Ion Reaction Cross 
Sections. Phys. Review C , vol. 41, no. 4, Apr. 1990, 

pp. 1610-1618. 

9. Ernst, David J.: Total Reaction Cross Sections and the 
Matter Density of Finite Nuclei. Phys. Review C , vol. 19, 
no. 3, Mar. 1979, pp. 896-904. 

10. DiGiacomo, N. J.; DeVries, R. M.; and Peng, J. C.: 
Microscopic Description of Nucleon-Nucleus Total Reac- 
tion Cross Sections. Phys. Review Lett., vol. 45, no. 7, 
Aug. 18, 1980, pp. 527-531. 

11. DeVries, R. M.; and Peng, J. C.: Nucleus-Nucleus Total 
Reaction Cross Sections. Phys. Review C, vol. 22, no. 3, 
Sept. 1980, pp. 1055-1064. 

12. Townsend, L. W.; and Wilson, J. W.: Comment on 
“Nucleus-Nucleus Total Reaction Cross Sections.” Phys. 
Review, ser. C, vol. 25, no. 3, Mar. 1982, pp. 1679-1681. 

13. Bradt, H. L.; and Peters, B.: The Heavy Nuclei of the 
Primary Cosmic Radiation. Phys. Review, second ser., 
vol. 77, no. 1, Jan. 1, 1950, pp. 54-70. 

14. Townsend, L. W.; Wilson, J. W.; and Bidasaria, H. B.: 
On the Geometric Nature of High-Energy Nucleus- 
Nucleus Reaction Cross Sections. Canadian J. Phys., 
vol. 60, no. 10, Oct. 1982, pp. 1514-1518. 

15. Towmsend, L. W.; and Wilson, J. W.: Energy-Dependent 
Parameterization of Heavy-Ion Absorption Cross Sec- 
tions. Radiat. Res., vol. 106, 1986, pp. 283-287. 


16. Townsend, L. W.; and Wilson, J. W.: Comment on 
“Trends of Total Reaction Cross Sections for Heavy Ion 
Collisions in the Intermediate Energy Range.” Phys. 
Review C , vol. 37, no. 2, Feb. 1988, pp. 892-893. 

17. Shen, Wen-Qing; Wang, Bing; Feng, Jun; Zhan, 
Wen-Long; Zhu, Yong-Tai; and Feng, En-Pu: Total Re- 
action Cross Section for Heavy-Ion Collisions and Its Re- 
lation to the Neutron Excess Degree of Freedom. Nucl. 
Phys., vol. A491, 1989, pp. 130-146. 

18. Vitturi, A.; and Zardi, F.: Modified Glauber Model for 
the Description of Elastic Scattering Between Heavy Ions. 
Phys. Review C , vol. 36, no. 4, Oct. 1987, pp. 1404-1407. 

19. Hussein, M. S.; Rego, R. A.; and Bertulani, C. A.: 
Microscopic Theory of the Total Reaction Cross Section 
and Application to Stable and Exotic Nuclei. Phys. Lett., 
vol. 201, no. 5, 1991, pp. 279-334. 

20. Menet, J. J. H.; Gross, E. E.; Malanify, J. J.; and Zucker, 
A.: Total- Reaction-Cross-Section Measurements for 30- 
60-MeV Protons and the Imaginary Optical Potential. 
Phys. Review C , vol. 4, no.- 4, Oct. 1971, pp. 1114-1129. 

21. Barashenkov, V. S.; Gudima, K. K.; and Toneev, V. D.: 
Cross Sections for Fast Particles and Atomic Nuclei. Prog. 
Phys., vol. 17, no. 10, 1969, pp. 683-725. 

22. Renberg, P. U.; Measday, D. F.; Pepin, M.; Schwaller, 
P.; Favier, B.; and Richard-Serre, C.: Reaction Cross 
Sections for Protons in the Energy Range 220-570 MeV. 
Nucl. Phys., vol. A183, no. 1, Mar. 20, 1972, pp. 81-104. 

23. Perrin, C.; Kox, S.; Longequeue, N.; Viano, J. B.; 
Buenerd, M.; Cherkaoui, R.; Cole, A. J.; Gamp, A.; 
Menet, J.; Ost, R.; Bertholet, R.; Guet, C.; and Pinston, 
J.: Direct Measurement of the 12 C + 12 C Reaction Cross 
Section Between 10 and 83 MeV/Nucleon. Phys. Review 
Lett., vol. 49, no. 26, Dec. 27, 1982, pp. 1905-1909. 

24. Buenerd, M.; Lounis, A.; Chauvin, J.; Lebrun, D.; 
Martin, P.; Duhamel, G.; Gondrand, J. C.; and 
De Saintignon, P.: Elastic and Inelastic Scattering of 
Carbon Ions at Intermediate Energies. Nucl. Phys., 
vol. A424, no. 2, Aug. 6, 1984, pp. 313-334. 

25. Kox, S.; Gamp, A.; Perrin, C.; Arvieux, J.; Bertholet, R.; 
Bruandet, J. F.; Buenerd, M.; Cherkaoui, R.; Cole, A. J.; 
El-Masri, Y.; Longequeue, N.; Menet, J.; Merchez, F.; and 
Viano, J. B.: Trends of Total Reaction Cross Sections for 
Heavy Ion Collisions in the Intermediate Energy Range. 
Phys. Review C, vol. 35, no. 5, May 1987, pp. 1678-1691. 

26. Dubar, L. V.; Eleukenov, D. Sh.; Slyusarenko, L. I.; and 
Yurkuts, N. P.: Parametrization of Total Cross Sections 
of Reactions in the Intermediate Energy Region. Sov. J. 
Nucl. Phys., vol. 49, no. 5, May 1989, pp. 771-773. 

27. De Vries, H.; De Jager, C. W.; and De Vries, C.: Nu- 
clear Charge-Density-Distribution Parameters From Elas- 
tic Electron Scattering. At. Data & Nucl. Data Tables, 
no. 36, 1987, pp. 495-536. 


13 


28. Bender, Carl M.; and Orszag, Steven A.: Advanced Math- 
ematical Methods for Scientists and Engineers . McGraw- 
Hill Boole Co., c. 19 78. 

29. Townsend, Lawrence W.: Harmonic Well Matter Den- 
sities and Pauli Correlation Effects in Heavy-Ton Colli- 
sions. NASA TP-2003, 1982. 

30. Bidasaria, Hari B.; and Townsend, Lawrence W.: Ana- 
lytic Optical Potentials for Nucleon- Nucleus and Nucleus- 


Nucleus Collisions Involving Light and Medium Nuclei . 
NASA TM-83224, 1982. 

31. Ringia, F. E.; Dobrowolski, T.; Gustafson, H. R.; Jones, 
L. W.; Longo, M. J.; Parker, E. F.; and Cork, Bruce: Dif- 
ferential Cross Sections for Small- Angle Neutron-Proton 
and Neutron-Nucleus Elastic Scattering at 4.8 GeV/c. 
Phys. Review Lett., vol. 28, no. 3, Jan. 17, 1972, 
pp. 185-188. 



4 6 

b, fm 


Figure 1. Absorption as a function of impact parameter for nuclear collisions. At half absorption, shaded areas 
5 l ~ 62- 



b, fm 

Figure 2. Absorption as a function of impact parameter for projectile-target systems of 4 He- 4 He and 16 0- 16 0 
at 2 GeV/nucleon. Terms R\ and R 2 denote geometric radii for systems. 
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Figure 3. Absorption as a function of impact parameter for three energy values: 50 and 250 MeV/nucleon and 
1 GeV/nucleon. Corresponding effective radii R\, /?2> and Rj are indicated. 



Figure 4. Variation of r 0 with atomic mass for harmonic oscillator density functions for nucleon-nucleus and | 

identical nucleus-nucleus collisions. . f 
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Figure 6. Variat 
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Figure 7. Variation of absorption cross section with energy for deuteron-nucleus collisions. 



i n L O I60-I60 




O Experiment 

Coulomb 

w/o Coulomb 


O 

I '0,0 



El, MeV/nucleon 

Figure 11. Absorption cross section as a function of energy for proton on carbon with and without Coulomb 
effects. 
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Figure 12. Absorption cross section as a function of energy for proton on aluminum with and without Coulomb 



Figure 13. Absorption cross section as a function of energy for proton on copper with and without Coulomb 
effects. 
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(a) Proton on lead. 



(b) Neutron on lead. 

Figure 14. Absorption cross section for proton or neutron on lead as a function of energy with and without 
Coulomb effects. Without Coulomb effects, curve corresponds to neutron projectiles. 
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Figure 17. Absorption cross section as a function of energy for oxygen on calcium with and without Coulomb 
effects. 
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Figure 18. Absorption cross section as a function of energy for oxygen on lead with and without Coulomb 
effects. 
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Figure 19. Absorption cross section as a function of energy for oxygen on silicon with and without Coulomb 
effects. 



Figure 20. Absorption cross section as a function of energy for calcium on calcium with and without Coulomb 
effects. 
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Figure 21. Absorption cross section as a function of energy for argon on bismuth with and without Coulomb 
effects. 



Figure 22. Absorption cross section as a function of [aJ ^ for protons on targets at 40 MeV/nucleon with 
and without Coulomb effects, where Aj is target mass number. 
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Figure 23. Absorption cross section as a function of ^A^/ 3 ^ for protons on targets at 60.8 MeV/nucleon with 
and without Coulomb effects, where Ap is target mass number. 



25 MeV/nucleon with and without Coulomb effects, where Ap is mass number of alpha particles and 
Ap is target mass number. 
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Figure 25. Absorption cross section as a function of ^Ay/ + Ap ^ for carbon on targets at 30 MeV/nucleon 
and 300 MeV / nucleon with and without Coulomb effects. 



2 

Figure 26. Absorption cross section as a function of + Ap^ for neon on targets at 30 MeV/nucleon 

with and without Coulomb effects. 
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Figure 27. Comparison of Bradt-Peters form with that for absorption cross section of proton on targets at 
60.8 MeV/nucleon. 



Figure 28. Comparison of Bradt-Peters form with that for alpha on targets at 25 MeV/nucleon. 
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Figure 29. Comparison of Bradt-Peters form with that for carbon on targets at 30 MeV/nucleon and 
300 MeV/nucleon. 



Figure 30. Comparison of Bradt-Peters form with that for neon on targets at 30 MeV/riucleon. 
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